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Abstract 

Consider the definition E of quasilocal energy stemming from the Hamilton- 
\ Jacobi method as applied to the canonical form of the gravitational action. We 

I \ examine E in the standard "small-sphere limit," first considered by Horowitz 

and Schmidt in their examination of Hawking's quasilocal mass. By the term 
small sphere we mean a cut S(r), level in an affine radius r, of the lightcone N p 
belonging to a generic spacetime point p. As a power series in r, we compute 
the energy E of the gravitational and matter fields on a spacelike hypersur- 
face S spanning S(r). Much of our analysis concerns conceptual and technical 
£f) ' issues associated with assigning the zero-point of the energy. For the small- 

sphere limit, we argue that the correct zero-point is obtained via a "lightcone 
■ reference," which stems from a certain isometric embedding of S(r) into a gen- 

uine lightcone of Minkowski spacetime. Choosing this zero-point, we find the 
following results: (i) in the presence of matter E = ^irr^lT^u^u^Wp + 0(r 4 ) 
"q I and (ii) in vacuo E = ■^r 5 [r^ u \ K u^u u u x u K ]\ p + 0(r 6 ). Here, is a unit, 

future-pointing, timelike vector in the tangent space at p (which defines the 
Sh . choice of affine radius); T^ u is the matter stress-energy-momentum tensor; 

T^,i,\ K is the Bel-Robinson gravitational super stress-energy-momentum ten- 
.£h ! sor; and \ p denotes "restriction to p." Hawking's quasilocal mass expression 

^ ' agrees with the results (i) and (ii) up to and including the first non-trivial 

C$ • order in the affine radius. The non-vacuum result (i) has the expected form 

based on the results of Newtonian potential theory. 
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INTRODUCTION 



Consider Einstein's non-covariant first-order action p|, the 4-integral of a "bulk" La- 
grangian which is quadratic in the Christoffel symbols and thus often called the "IT ac- 
tion." Starting with the Einstein action, one applies standard techniques associated with 
Noether's theorem in order to derive, among other things, an energy definition in general 
relativity: namely, the 2-integral of an Einstein superpotential over some generic 2-surface 
S in spacetime M.Q The Einstein energy is well known to be ambiguously defined because 
it depends on the choice of background coordinates. Nevertheless, one may use the Einstein 
construction to define sensible notions of total gravitational energy. Indeed, consider the 
scenario of asymptotic flatness, say, towards future null infinity J + . In this case, S tends 
to a round, infinite-radius, 2-sphere cut of J + , and the (now suitably unique) choice of 
asymptotically Cartesian coordinates ensures that the Einstein energy agrees with the ac- 
cepted Trautman-Bondi-Sachs (tbs) notion of total energy. || However, were we to offer the 
Einstein definition as the energy contained within some quasilocal (that is, finite) 2-surface 
S, we would still be confronted with the task of choosing a physically meaningful set of 
background coordinates. The only natural choice would be coordinates which are partially 
adapted to the embedding of S C M. However, such a choice wrecks the agreement between 
the Einstein and TBS energies as S tends towards J + . In fact, choosing such S'-adapted co- 
ordinates, one finds that the Einstein energy blows up in the said limit. Similar statements 
can be made regarding other approaches for defining energy, which trade coordinate (that is, 
holonomic-frame) ambiguity for ambiguity of a different stripe, e. g. tetrad (or rigid) frame, 
spin frame, or auxiliary vector (or spinor) fields. The traditional party line regarding these 
issues is the following: there is no over-arching rule, applicable for all quasilocal 2-surfaces, 
for selecting a (suitably unique) background frame; whence gravitational quasilocal energy 
is not well-defined.Q To what extent does the stubborn presence of frame ambiguity in the 
quasilocal context point to a gap in our understanding of gravitational energy? To address 
this question and sharpen our thoughts on these issues, let us consider a covariant version 
of the Einstein construction. 

Employing a straightforward field-theoretic generalization of the Hamilton- Jacobi (hj) 
method ||, one may derive from a covariant action functional an expression for canonical 
quasilocal energy (qle) in general relativity. We call this definition of QLE canonical, 
because, owing to its intimate connection with HJ theory, this QLE is also the on-shell value 
of the gravitational Hamiltonian for the choice of unit lapse function and vanishing shift 
vector at the system boundary. || We also note that the canonical QLE is the thermodynamic 
internal energy in a thermodynamical description of a (relativistic) self-gravitating system. 



1 See the excellent review article by Goldberg in Ref. for the details of this analysis. 

2 A serious contender for such an over-arching rule has been given by Dougan and Mason, who 
use certain "(anti-)holomorphic" spinor fields in order to define a four-dimensional space of "quasi- 
translations" pointing on and off an essentially generic 2-surface. || Szabados has shown that the 
Dougan-Mason proposal provides a tidy characterization of pp-wave geometries. 
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The analysis that leads to the canonical QLE runs along a somewhat different line than 
the one followed in a Noether-type analysis, but it also leads to a concept of energy which 
is not unique. Indeed, as is always the case with energy, the canonical QLE is defined only 
up to the choice of a zero-point. The zero-point ambiguity may be traced to a freedom 
present in the action principle. Namely, one may always add to the action any functional of 
the fixed boundary data without affecting the variational principle. As with the situation 
above, if the goal is to obtain agreement with the accepted notions of gravitational energy 
at spacelike or null infinity, then there is a suitably unique choice of energy zero-point ||[8| , 
whereas at the quasilocal level there seems to be no preferred choice.^ While at first sight 
this seems no better or worse than the situation encountered above, notice that now the 
ambiguity in the energy has a physical interpretation, and, moreover, is a field-theoretic 
generalization of the standard ambiguity present in the HJ definition of energy in ordinary 
mechanics. We may now restate the emphasized portion of the party line above as follows: 
there is no over-arching rule, applicable for all quasilocal 2-surfaces, for selecting a (suitably 
unique) zero-point. Taking this statement at face value, we claim that it is the physicist's 
job to select an appropriate energy zero-point, guided by the principle that the selection 
should be appropriate for the physics of the problem at hand. We would like to point out 
that this is a common enough state of affairs in general relativity, a many-faceted theory 
known for its wealth of possible boundary conditions. Indeed, by way of analogy consider 
the search for solutions of the Einstein field equations. In practice, relativists certainly do 
not attempt to find the general solution, rather they attempt to find solutions given some 
additional physical input (boundary conditions, symmetries, etc.). In practice, the same 
such additional input is needed to associate a meaningful QLE with a particular quasilocal 
2-surface. 

Bearing these points in mind, we recall the form of the canonical QLE: 



where we adopt geometrical units (in which both Newton's constant and the speed of light 
are set to unity), S is a closed 2-surface, dS is the proper area element on S, and k is the 
mean curvature of S as embedded in a spanning 3-surface S. It is important to realize that 
this E, while obtained as a proper surface integral over S, is the energy of the gravitational 
and matter fields which live on S, that is to say, E is a functional of the initial data of E. 
This concept of energy is rooted in the 3 + 1 view of spacetime geometry, and for a fixed S it is 
slightly sensitive to the choice of spanning S. [More precisely, E depends on the equivalence 
class of spanning 3-surfaces determined by a unit timelike vector on (and pointing orthogonal 
to) S.} This sensitivity is quite analogous to the observer dependence of energy in special 
relativity, and a priori we expect its presence. || In general, the term k\ Ict = k\ rcf (a a b) 



3 Other than the choice of a vanishing zero-point. This choice corresponds to the aforementioned 
choice of S"-adapted coordinates for the Einstein definition. Like before, such a choice, leading to 
an infinite energy, wrecks the agreement between the canonical quasilocal and TBS energies in the 
large-sphere null limit. 




(0.1) 
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represents an arbitrary (local) function of the intrinsic metric a a b of S and corresponds to 
the freedom to assign the QLE zero-point. Notice that this freedom corresponds to a proper 
surface integral of what is effectively a free function of two variables. This freedom, stemming 
from the field-theoretic character of gravity, is rather more subtle than the freedom in simple 
mechanics of simply adding a constant to the energy. For our analysis here we only consider 
two energy zero-points, one determined by lightcone reference and the other by Euclidean 
reference. We make these concepts precise below. 

In this paper we examine the canonical QLE in the standard "small-sphere limit," first 
considered by Horowitz and Schmidt in their classic examination of Hawking's quasilocal 
mass.Q By the term small sphere we mean a cut S(r), level in an affine radius r, of the 
lightcone N p belonging to a generic spacetime point p. As a power series in r, we compute the 
energy E of the gravitational and matter fields on a spacelike hypersurface E spanning S(r). 
Much of our analysis concerns conceptual and technical issues associated with assigning the 
zero-point of the energy, and, therefore, particularly elucidates the points raised in the first 
two paragraphs above. For the small-sphere limit, we argue that the correct zero-point is 
obtained via the aforementioned lightcone reference, which stems from a certain isometric 
embedding of S(r) into a genuine lightcone of Minkowski spacetime and amounts to fixing 

fc| ref = -y / ±ft{l + |(ft)- 1 Alogft + [1 + 1(7^-^ log ft]" 1 } (0.2) 

in the above QLE definition. Here A is the Laplacian operator and TZ is twice the Gaussian 
curvature of S(r); therefore, as advertised, this choice for fc| ref depends solely on the intrinsic 
geometry of S(r). Notice that, due to the presence of the square root in this expression, one 
expects this choice for k\ re{ to be valid only for (topologically spherical) 2-surfaces possessing 
everywhere positive Gaussian curvature (as is the case both for the small spheres we study 
here). Choosing the proper surface integral of this choice for (87r) _1 /c| ref as the energy zero- 
point, we find the following results for small spheres: in the presence of matter 

E = l7cr 3 [T^uV]\ p + 0(r 4 ) (0.3) 

and in vacuo 

E = ±r 5 [T^ XK u»u u u x u K ]\ p + 0(r 6 ) (0.4) 

Here, is a unit, future-pointing, timelike vector in the tangent space at p (which defines 
the choice of affine radius); T^ u is the matter stress-energy-momentum tensor; T^ak is the 
Bel-Robinson gravitational super stress-energy-momentum tensor;^ and \ p denotes "restric- 
tion to p." It is interesting to note that, when integrated, the Bel-Robinson "energy" in ( p.4| ) 



4 This limit is also considered in Ref. [11] by Kelly, Tod, and Woodhouse for Penrose's kinematic 
twistor and associated quasilocal mass, and in Ref. [12| by Dougan for the Dougan-Mason quasilocal 
four-momentum. 



5 Our vacuum-case Bel- Robinson tensor is the following: 
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has been proven to be very useful already in the sense of a mathematical "energy" in the 
study of existence of solutions of hyperbolic equations in Einsteinian gravity. [13-[L6j] It is 
noteworthy that it shows up also in the physical limit given in ( p.4| ). Although the full phys- 
ical significance of the physical limit is not known to us, we think that these mathematical 
and physical properties go beyond mere coincidences. For both the vacuum and non- vacuum 
cases, Hawking's quasilocal mass expression agrees with the canonical QLE results (|0T3| ) and 
( p.4| ) up to and including the first non-trivial order in the affine radius. We find this result 
rather striking in light of the fact that the Hawking mass has no apparent connection with 
the gravitational action or Hamiltonian. We show that the non-vacuum result ( P-ci| ) has the 
expected form based on the results of Newtonian potential theory. 

Compare the small-sphere limit considered here with the large-sphere limit when S(r) 
tends to a cut of J + , in which case we know that the choice of Euclidean reference yields 
agreement between E and the TBS energy. || In both limiting cases, the 2-surface S(r) of 
interest is a cut, level in an affine radius r, of an outgoing congruence of null geodesies. There 
is, however, a crucial distinction to be made. In the small-sphere case, S(r) arises as the 
cut of a genuine lightcone, while in the large-sphere case this is generally not true. We find 
it remarkable that this distinction can be mirrored in the choice of zero-points. To grasp 
this point, consider first the Euclidean reference for either limit. This reference involves 
an isometric embedding of S(r) into a flat Euclidean 3-space E 3 . Now, in both limiting 
scenarios S(r) is, in general, slightly distorted from perfect roundness. Therefore, with E 3 
viewed in turn as an inertial slice of Minkowski space M, the enveloped S(r) C E 3 C M 
cannot be the cut of a genuine lightcone of M. (Technically, in this case the outward null 
congruence associated with S(r) is not shear-free, but the lightcones of M are shear-free.) 
Therefore, in the general large-sphere scenario, neither the S(r) embedding into the physical 
spacetime M nor its embedding into the reference spacetime M corresponds to a lightcone 
embedding [that is, in neither case is S(r) the cut of a lightcone]. Of course, for the small- 
sphere limit we may employ either the Euclidean reference or the lightcone reference. Of 
these two choices, the lightcone employs flat spacetime to put the reference space on an equal 
footing with the small sphere construction in the physical spacetime. To define the lightcone 
reference, we first isometrically embed S(r) into the lightcone of a point q G M, and then 
select a certain 3-surface T C M spanning S(r) [so that S(r) = Hf]N q }. The details of 
this construction, eventually leading to the expression (p.2|), are found in Subsection 2. A. 
Now, having defined the lightcone reference (tailored to the small sphere limit) and found 
the resulting closed-form expression ( p.2[ ), we may now invert the question, asking whether 
or not a lightcone- reference fc| rof defines a QLE (p.l|) possessing a correct large-sphere limit 
towards J + . We discuss this issue in Appendix A. 

The organization of this paper is as follows. In Section 1 we lay the foundations for our 



where C^pXa is the Weyl tensor, *C IMp \ (7 := ^e^Q/jC^Aff is the left-dual of the Weyl tensor, and 
± Cpp\a '■= \{C 'npXaT^itfipapC ^ Ao-) is the self-dual (+) [anti-self-dual (— )] part of the Weyl tensor. 
Further curvature conventions are discussed in the appendix. 
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examination of the small-sphere limit. We describe the geometry of the limit in Subsection 
l.A, fix some general conventions in Subsection l.B, and make some general observations 
concerning the embedding of a 2-surface in Minkowski spacetime in Subsection l.C No choice 
of energy zero-point is made in Section 1, although the results of Subsection l.C are used 
in the subsequent sections to construct zero-points. In Section 2 we study the small-sphere 
limit, subject to the choice of lightcone reference, and derive the results (|0.3j ) and ( p.4|) . 
In Section 3 we discuss the relationships between the main results (p.3|), ( p.4[ ) and results 
from Newtonian potential theory. Appendix A contains discussions of the small-sphere limit 
subject to the choice of Euclidean reference, and the large-sphere limit towards J + subject 
to the choice of lightcone reference. Throughout our analysis, we use the Newman- Penrose 



(np) formalism |T7|)|H1I? w hh which we assume the reader is familiar. In Appendix B we 
collect various conventions and results associated with the NP formalism which are used in 
the main parts of the paper. 



I. PRELIMINARIES 
A. Geometry of the limit construction 

Choose a generic spacetime point p e M, as well as a unit, future-pointing, timelike vector 
m m which lies in T p (M), the tangent space at p. We may think of as the instantaneous four- 
velocity of an Eulerian observer at p. Let N p C M represent the future lightcone generated 
by the null geodesies emanating from p. Label the generators of N p by coordinates (9,(j)), 
or equivalently by (CjC)j w here £ := e 1 ^ cot (9/ 2) is the stereographic coordinate. Further, 
choose the afline parameter r along the generators of the lightcone which at the point p 
satisfies the following conditions: r = and the null tangent Z M := (<9/<9r) M to the aflinely 
parameterized generators obeys Fu^ = —1. By the term small sphere we shall mean a 2- 
surface S(r) C N p level in the coordinate r. Provided that we restrict our attention to small 
enough values of r, we need not be troubled by conjugate points and each 2-surface S(r) 
will be only slightly distorted from perfect roundness. On our lightcone N p we construct 
a null tetrad n M , m M , m M } as follows. We take rF as the inward null normal to each 
S(r) (normalized so that Fn^ = —1), and assume that the complex space leg m M points 
everywhere tangent to each S(r). Further, via enforcement of the condition m A / M V M mA = 0, 
we remove the freedom to perform r-dependent rotations of the complex dyad. Together 
with the geodetic property of Z M , this implies that the spin coefficient e = [cf. Eq. ([B"T|a)l. 

We also consider a standard pseudo-orthonormal tetrad {u^ =: e_i_ M , e^, e y M , e z M } at the 
point p, in terms of which we have the following expansions at p: 

F = m m + sin#cos0e x M + sin#sin0e y M + costfe/ (1.1a) 
n» = | [u 11 - (sin 9 cos 4> e/ + sin 9 sin 4> e/ + cos 9 e/)] . ([TT|b) 

Capital Latin letters, e. g. A, B,C,- ■ ■, denote pseudo-orthonormal tetrad indices and run 
over the values (_L,x, y, z). The expansions F = I a ca^ and ra M = n A eA^ in Eq. ( |1 . 1| ) show 
that the components l A and n A are essentially the first four spherical harmonics. Standard 
orthogonality properties of the spherical harmonics then yield the following identities: 
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d(0) l A l B = lu A u B + \g AB (1.2a) 

d(Q.) l A n B = \u A u B - \g AB , QT2|b) 

J d(Q) l A l B fl D = \ [16u a u b u c u d + 12u( A u B g CD > + g {AB g CD) ] , Q) 

which prove quite useful for reducing many of the integral expressions encountered below. 
Here we make use of the convenient notation d{Q) := (47r) _1 dfi, where dfl is the area element 
of a unit-radius round sphere and the integrations in Eq. (|1.2j ) are over the unit-radius round 
sphere. 



B. Physical and reference energy surface densities 

Let us now define the physical energy surface density (87r) _1 A;, whose proper surface 
integral is the total unreferenced quasilocal energy. Our construction requires that we fix 
a 3-dimensional hypersurface £ spanning S(r), or, more precisely, an equivalence class of 
spanning 3-surfaces determined by the choice of a unit, future-pointing, timelike vector w M 
on (and pointing orthogonal to) S(r). We choose 

Ip + n", (1.3) 

which at p agrees with the four-velocity introduced in the last paragraph. In terms of the 
convergences p and —/j, of the null normals defined in the appendix Eqs. (|BT| i,k), the mean 
curvature of S(r) as embedded in £ is given by 

k = 2p + p. (1.4) 

We shall write 

E\ ph * := (87T)- 1 I dSk (1.5) 

J.S(r) 

for the unreferenced QLE associated with the physical slice S in spacetime M. Likewise, we 
introduce the reference energy surface density (87r) _1 fc| rof , and with it define the reference 
contribution to the QLE, 

£| rcf : = (87T)- 1 I dSk[ ct . (1.6) 

Js(r) 

As mentioned, E\ Ici is a functional solely of the intrinsic geometry of S(r), although as yet 
we have made no definite choice for this functional. In Section 2 we choose the specific 
functional stemming from the lightcone reference, while in Section 3 we choose the one 
stemming from the Euclidean reference. The difference 

E = E\ phy - E\ Ici (1.7) 

is the total referenced QLE ( |0.1| ) . 
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C. The geometry of 2-surfaces in Minkowski spacetime 

This subsection collects a few basic results concerning the reference embedding S(r) C M. 
We shall use these results later when constructing a particular reference energy surface den- 
sity. For notational convenience here and in what follows, we often use a sans serif notation 
for objects associated with Minkowski spacetime M, and we may write k in place of k\ re{ . 
We restrict our attention to Minkowski- spacetime references. That is to say, the reference 
energy surface density (87r) _1 k is determined via an auxiliary isometric embedding of the 
2-surface S(r) into a 3-dimensional hypersurface Z which is itself contained in Minkowski 
spacetime M. Physically, this means that we assign the zero value of the energy to that 
portion of the slice Z contained within S(r). 

At the end of the day, our expressions for k depend only on the S(r) 2-metric o~ a b 
(with a, b, c, ■ ■ ■ as S(r) indices), and, therefore, there is technically no need to consider the 
reference spacetime. Nevertheless, in order to motivate and derive our choices, we begin 
with the chain of inclusions S(r) C Z C M and an overall spacetime point of view associated 
with it. Let us collect a few definitions. Take as the unit future-pointing normal of Z 
as embedded in M. Let k a b be the extrinsic curvature tensor associated with the spacelike, 
outward-pointing, unit normal of S(r) as embedded in Z. Denote by the extrinsic 
curvature tensor (with i,j, k, ■ ■ ■ as Z indices) of Z associated with u M . Projection into S(r) 
of all of the free indices of K^- defines an extrinsic curvature tensor l a & on S(r) associated 
with u M . Tangential-normal projection with respect to S(r) of Kjj defines a covector on 
S(r). If S(r) is to ctriSG clS db 2-surface in M, then along with its intrinsic metric a a b the 
triple {k a b, \ a b, A&} of extrinsic data must obey certain constraints. These are integrability 
criteria relating the intrinsic and extrinsic data of S(r) to the (vanishing) components of the 
M Riemann tensor. Among the constraints for S(r) C M are the following:^ 

(k) 2 - k ab k ab - (I) 2 + \ ab \ ab -TZ = (1.8a) 
Kb - k\ a + Ad - A a l a6 = . dTUb) 

Here JZ is the scalar curvature of S(r), and the colon denotes covariant differentiation in- 
trinsic to S(r). We shall not have need to consider the other embedding constraints. 

II. SMALL-SPHERE LIMIT 

Throughout this section, the term A;| rcf stands for the explicit expression ( p.2| ) given in 
the introduction, and the functional ( |1.6|) is fixed accordingly. 



6 As our use of power series in r is widespread in this paper, we use parenthesis when a variable, 
say 1Z, is raised to a power. Hence, (TZ) 2 = 1Z1Z and (TZ)^ 1 = 1/TZ, while 1Z 2 always denotes the 
0(r 2 ) coefficient in the expansion 1Z = ■ ■ ■ + r 2 lZ 2 + • • - . The only exception to this rule will be the 
radius r itself. As there is no possibilty for confusion, we use, for example, r~ 1 to mean 1/r. 
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A. Lightcone reference 



To construct the lightcone reference (see the figure at the end of this paper), consider the 
notations of Subsection l.C and assume that S(r) is a cut of a genuine lightcone N g belonging 
to a point q G M. In this case, the geodesic congruence N g is sheer-free, which means that 
the complex shear <; := k mm + \ mm of N g vanishes. Here k mm := k ab m a m b and \ mm := \ a bm a m b 
are complex components with respect to the S(r) dyad {m a , m a }, respectively capturing the 
trace-free pieces of k ab and \ ab . The vanishing of q thus implies that the trace-free piece of 
k ab equals minus the trace-free piece of \ ab ; whence Eq. (pLffit) becomes 



(k) 2 - (I) 2 -211 = 0. (2.1) 

We have no need to consider Eq. ( |1.8| b) in this subsection.[] We define the convergence 
g := ~ (k — I) of the outward null normal to S(r) C N g , and with it rewrite Eq. (|2.1|) as 

(k) 2 -(k-2g) 2 -211 = 0. (2.2) 

Now, we are going view g = g((, £; r) as some specified function. Specification of g (indi- 
rectly) fixes a slice E (or, more precisely, an equivalence class of slices) spanning S(r), which 
in turn defines k. With the last equation and some algebra, we find 

k=\{g)- l K + g (2.3) 

as the expression for k in terms of the function g which is not yet specified. Now, plugging 
both the radial expansion (|B21|) for 1Z and some radial expansion for g into Eq. Q2.3|) , we 



obtain a radial expansion for k, concerning which we make the following crucial observations. 
First (relevant for the non- vacuum case), a reference convergence of the form g = — r _1 +0(r) 
determines an expansion for k up to and including 0(r 2 ) (which is actually higher than 
the order needed to get the non- vacuum limit). Second (relevant for the vacuum case), a 
reference convergence of the form g = — r _1 + 0(r 2 ) determines an expansion for k up to 
and including 0(r 3 ). Since we wish our final expression for k to serve as a proper reference 
term, we demand that g = g((, (; r) be built purely from the 2-metric a ab of S(r) (so that k 
is). This restriction alone hardly fixes the choice of g (and k). However, for both scenarios 
of interest the geometry of the embedding of S(r) into the physical spacetime M suggests a 
natural choice. 

Let us write down our choice for g, verify that it leads to Eq. (p.2|) , and finally discuss 
why it is physically meaningful. We pick 



g = -J\K [1 + K^r^cllogft] , (2.4) 



7 For the lightcone reference we construct in this subsection, Eq. (|1.8| b) and the other embedding 
constraints not appearing above would be differential equations determining the remaining S(r) 
extrinsic data from a ab and our choices for I and k. 
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where 9 is the full "eth" operator on S(r) [cf. Eqs. ( |B16|) and (P14j) 1 and on weight-zero 



scalars the S(r) Laplacian is A = 259. With this choice for g (and hence this choice for E) 
we obtain from Eq. (|2.3|) the same closed-form expression, 



!ft{i + !(ftr 1 6giogft+ [l + K^-^aiog^]- 1 } , (2.5) 

as given in Eq. ( p.2| ) for 811 times the reference energy surface density. Now, let us argue that 
Eq. fl2.5|) defines a physically sensible reference surface density for the small-sphere limit. 
Using Eqs. ( |B16| ) and ( [B14| ) for the <3 operator and the radial expansion ( |B21| ) for 7Z, we 



find (whether in vacuo or not) that (|2.4j) satisfies 

Q = -r- 1 - \r{TZ° + \%$ Q K Q ) + 0(r 2 ) , (2.6) 

where 9 is the unit-radius round-sphere "eth" operator. Hence, we have at least g = 
— r _1 + O(r) for the reference convergence. To start with, consider the non- vacuum case and 
glance at the expansion (|R4]a) for the convergence p associated with the physical embedding 
S(r) C N p . We see that through 0(r°) our choice for g agrees with the physical p. Now turn 
to the vacuum case and notice that the expansion ((B5|a) for the convergence p associated 
with the physical embedding obeys p = — r _1 + 0(r 3 ). Combining Eq. (|2.6|) with the vacuum 
identity^ 

g 6 ^° = -3TZ , (2.7) 

we have g = — r~ l + 0(r 2 ); and, therefore, now through Ofr 1 ), our choice for g agrees with 
the physical p. To sum up, we can state that, whether in vacuo or not, our choice ( |2.4| ) 
for g determines an embedding of S(r) into N 9 C M which would seem closely related both 
intrinsically and extrinsically to the physical embedding of S(r) into N p C M. Let us now 
put aside the issue of motivation, and simply expand our choice (|2.5|) for k in powers of r. 
Since g = — r _1 + O(r) in non- vacuum, from Eq. ( |2.3| ) we get 

k\ rct =: k = -Ir- 1 - \rK° + 0(r 2 ) (2.8) 

for the non- vacuum radial expansion of Eq. (|2.5|) . Likewise, since g = — r _1 + 0(r 2 ) in vacuo, 
from Eq. ( ^.3| ) we now get 

=: k = -2r~ l - \rV? - \r 2 U l - \r z V? + 0(r 4 ) (2.9) 

for the vacuum radial expansion of Eq. ( |2.5| ) . 



8 One proves the identity as follows. First, as we do in the appendix, calculate the vacuum-case 
0(r°) piece of the S(r) curvature scalar. One finds TZ° = -4(<J^ + ^) [cf. Eqs.(jB25|a) and (p24|d)]. 
But, as shown in the final part of the appendix, the radial Bianchi identities [cf. Eq. QB24 d,e,f,g)] 
imply that SoSo^^ 



-Sty® which establishes the result. 
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B. Non- vacuum limit 



Turning our attention to the non- vacuum scenario, we begin by using the appendix 
Eqs. (P^a,f) for the convergences —\l and p of the null normals to determine the following 
expansion for 87T times the physical energy surface density: [cf. Eq. (|1.4j )1 

k = ~2r- 1 + r + ^° + 2A° + |$°J + 0(r 2 ) , (2.10) 

Here, ^/^ ^ii> and ^° are r — * limits of standard curvature terms from the NP formalism. 
Respectively, they are proportional to a certain null-tetrad component of the Weyl tensor 
C^uXk at p, a certain null-tetrad component of the Ricci tensor at p, and the Ricci scalar 
R at p [cf. Eqs. ( |B3| 3,f,k)1. Because we have adapted our null tetrad to the lightcone N p , 
the components ^ and ^11 are m fact angle-dependent, i. e. they are functions of (C C)j 
however, as the limit of a scalar function, A is not angle dependent. Next, we substitute 
both Eq. ( p. 10] ) and the appendix Eq. ( |B18|) for the surface element dS of S(r) into the 
basic expression (|1.5|) for E\ pby , and expand out the result in order to obtain 



E \^y = _ r + r 3 J d ( ) [i(^o + t 0) + A o + i $ q + i$oj + 0(r4) (2 n) 

Like before, <3>q , an angle-dependent term, is proportional to a certain null-tetrad component 
(|B"3|g) of the Ricci tensor at p. We can perform the integrations remaining in Eq. Q2.11| ). 
The spherical average of the real part of ^S] vanishes identically. Indeed, the r — > limit 
of Eq. (|B3|c) shows that ^ + ^2 = (l A n B l c "n L ) C 'abcd)\ p , and by symmetry the average of 
l A n B l c n D must yield terms either proportional to u A u B u c u D or containing g AB . Therefore, 
that the average in question vanishes follows from the index symmetries and trace-free 
character of the Weyl tensor. For the next integration we can straightaway write J d(O) A = 
jgR\ p - To evaluate the unit-sphere averages of the other terms in the integrand, we must use 
the identities (Qa,b). The r -> of Eq. (g^g) is <£>[J = \[l A l B R AB \\ p , which shows that the 
average of $q is readily obtained with Eq. ( |1.2| a). By rewriting the r — > limit of Eq. (|B3| k) 
as = [\I a ti b Rab + ^R]\p, w e likewise obtain the average of with Eq. (|0|b). Adding 
together the individual results for these integrations, we find 

E\** = -r + ^r 3 [bR^vr + 2R)\ p + 0(r 4 ) (2.12) 

as the final limiting expression for the unreferenced qle. 

To obtain an expansion analogous to Eq. ( j2.12j ) for E\ TcC , we must first explicitly compute 
the radial expansion for 8ir times the reference energy density. Putting together Eqs. ( p.8|) 
and (|E2^ ), we get 

k\«* = -27-- 1 + |r [25 *° + 29 *° - S $° - 8 <l? - $° ] + 0(r 2 ) . (2.13) 



Now, substitute Eqs. (|2~l3| ) and ( [B18D into Eq. Q), do some algebra, and perform a few 



simple integrations^] in order to reach 



9 Here and in what follows, "simple integrations" refer to the following version (and its complex 
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E\ 



-- -r + ±r 3 J d(n) $° 00 + 0(r 4 ) . (2.14) 

Like before, we use Eq. ( |1.2| a) to compute the remaining integral, thereby finding the desired 
result: 

= - r + M r3 [*R^u v + R]\ p + 0(r 4 ) , (2.15) 

which should be compared with Eq. ( |2.12j ). 

We may now easily obtain the Introduction's result (p.3|) for the non-vacuum limit. 
Indeed, combination of Eqs. ( |2.12D and ( |2.15| ) yields the following result for the total QLE 
(IT 



£ = ^[G^Y]| p + 0(r 4 ). (2.16) 

Notably, this result for E is valid whether or not the field equations hold, i. e. it is a 
geometric identity. However, with Einstein's equations G^ u = StcT^ we immediately arrive 
at the Introduction's result ( |0.3| ). 



C. Vacuum limit 

The derivation of the vacuum limit proceeds along the same lines as those just con- 
sidered for the non- vacuum limit. With the appendix Eqs. (p5|a,f) for the vacuum-case 
convergences — /i and p, we compute the radial expansion for the physical energy surface 
density, [cf. Eq. flOP ] 

k = -2r~ l + \t [S *S + S *?] + \r 2 [S *{ + S *J] 

+r 3 [&>*°o*o - io^o (*8*? + 4*?) - ±5 + 4*?) - \n*] + 0(r 4 ) . (2.17) 

Next, we substitute this expansion along with the expansion (|B19| ) for the vacuum-case 
surface element into Eq. ( |1.5| ), do some algebra, and perform a few simple integrations. 
These steps lead to 

Er = -r + r 5 J d<0> [^*°*° + \ - V})\ + 0(r 6 ) . (2.18) 

A glance at the explicit expression ( |B23| c) for the coefficient TZ 2 [of the 0(r 2 ) term in the 
expansion ( |B21| ) for TZ] shows that the unit-sphere average of the term within the parenthesis 
above vanishes. To evaluate the final unit-sphere average, we first note that the r — > limit 
of the square of Eq. (|B3]a) is in fact *°*° = [\I a I b I c I d T ABC d\\p- This may be checked by 



conjugate) of the divergence theorem. Suppose that sw(/) = — 1 (with sw denoting spin weight). 
Then with 8 representing the "eth" operator on the unit-radius round sphere, we have J d(O) S / = 
0. 
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an explicit computation using our definition ( |0.5| ) of the Bel-Robinson tensor. Evidently 
then, we may use Eq. flL2|c) to find 



E\ phy = -r+ ^r 5 [Tr, pa vru"u p u a \\ p + 0(r 6 ) (2.19) 

as the desired limit expression for E\ phy in vacuo. 

Turning now to the calculation of E\ Icf , we first put together Eqs. ( |2.9|) and ( B23|) , in 
order to get the following explicit expansion: 

fc |»f = _ 2r -i + | r p o ^o + g^O] 

+^r 2 [a *l + 3 *i] + |r 3 - TZ 2 ) - ^*°*°] + 0(r 4 ) . (2.20) 



We substitute this expansion as well as the expansion ( pi9|) into Eq. ( |1.6|) , and again do 
some algebra and a few simple integrations, thereby obtaining 

E\^ = -r + r 5 j d(Q) [^jjfcg + \ (4*8*8 ~^)]+ 0(r 6 ) . (2.21) 
Finally, calculations identical to those just performed for E\ pby establish that 

E \" = ~ r + 2l5 r5 [W^^llp + °( r6 ) > ( 2 - 22 ) 

and the difference of ( |2.19| ) and ( |2.22|) immediately gives the Introduction's result (|0.4|). 



III. NEWTONIAN POTENTIAL THEORY 

It is interesting to compare our main results with analogous results from Newtonian po- 
tential theory.0 The Newtonian interpretation for the Introduction's result (p.3|) is straight- 
forward. Consider a pressureless ball of fluid of radius r and constant (volume) mass density 
e. The total Newtonian mass for the ball is M = |7rr 3 e. If we identify this radius with 
the affine radius that appears in the result (|0.3| ), then a further identification between the 
Newtonian mass M and the quasilocal energy E yields the correspondence e = [T^u^u"} \ p . 
This is a solid result, fully in accord with the non-existence of "pure gravitational field en- 
ergy" at a point, since only matter energy contributes in this limit. The Newtonian analog 
of the result ( p.4j) is not clear. Nevertheless, we find the following result rather interesting. 
Consider again the ball of fluid, but now let us compute the Newtonian gravitational field 
energy within the ball. Recall that the energy density of the Newtonian gravitational field is 
given by — V$- V$/(8tt), where V is the flat-space gradient operator and $ is the Newtonian 
potential. Inserting the appropriate expression for $ into the energy density and integrating 
over the interior of the ball, we obtain 

E N = -^r 5 (47re) 2 (3.1) 



10 A comparison between the canonical QLE with Euclidean subtraction and the Newtonian grav- 
itational energy is given in Ref. 0. 
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for the gravitational potential energy inside the ball. Comparison of this result with the 
small-sphere result ( p.4|) , we see that both expressions depend on the fifth power of radius and 
both contain a numerical factor of ^. We do not know at this time if the close resemblance 
between these results has any physical significance. 
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APPENDIX A: EUCLIDEAN REFERENCE AND LARGE-SPHERE LIMIT 

One can also carry out our analysis of the small-sphere limit for the choice of Euclidean 
reference. To define the Euclidean reference, again consider the notations of Subsection I.C, 
but now assume that Z is a flat inertial hyperplane E 3 in M. In this case, \ ab = = A b and 
the constraints given in Eqs. (|1.8| a,b) are simply 

(k) 2 - k afe k a6 -TZ = (Ala) 
k: 6 -k° 6:a = 0. (0b) 

To analyze these equations, it proves useful to split k ab into its trace and trace-free pieces. 
This is readily achieved by working with the components of k ab with respect to the null dyad 
{m a ,rh a }. Indeed, 2k m?Ti := 2m a m b k ab and k mm := m a m b k ab respectively capture the trace 
and trace- free pieces of k ab , and in terms of these quantities the Eqs. (|Al| a,b) become 

[^m?fi] kmj^kjf^Tj 2^ ^ (A2a) 

<3k m m — (3k m m = . ( |A2| b) 

For the case of Euclidean reference, we are unable to obtain a closed-form expression for k 
(neither in the non- vacuum nor vacuum cases). Nevertheless, it is evident from Eqs. (|A"l|a,b) 
that the full extrinsic curvature tensor k ab (and hence its trace piece) is determined solely 
by the intrinsic metric on S(r). One expects that Eqs. (|A2|a,b) may be solved for k ab , 
provided S(r) is only slightly distorted from perfect roundness. For small enough values of 
r, the expansion ( |B21| ) for TZ assures us that this is indeed the case. Moreover, the solution 
k ab should be unique up to Euclidean translations and rotations. Lacking a closed-form 
expression for k, we have obtained radial expansions in both the non-vacuum and vacuum 
scenarios. For both scenarios, we obtain such expansions by first plugging into Eqs. (|A"2|a,b) 
the Ansatze 



14 



kmm — r + rk mfh + r k m}fl + r k mm + 0(r ) (A3a) 
k m m = rk mm + r k mm + r 3 k^ m + 0(r ) , ( |A3| b) 

along with appropriate (vacuum or non-vacuum, as the case may be) expansions for both <3 
and TZ. We have performed this calculation, and find that the choice of Euclidean reference 
also establishes ( |0.3| ) in the non-vacuum case. However, with the choice of Euclidean refer- 
ence, the vacuum small-sphere limit of the QLE is not related directly to the Bel-Robinson 
tensor. 

Our closing comments concern the large-sphere limit, in which case S(r) tends to a 
round infinite-radius 2-sphere cut of J + . Thus, we now consider a spacetime M which is 
asymptotically flat towards J + and a corresponding system of Bondi coordinates (w, r, (,(,)■ 
Here w is a retarded time coordinate, and r is an affine radius similar to before. Now S(r) 
arises as a cut, level in r, of an outgoing null hypersurface M (in general not a genuine 
lightcone), and twice the Gaussian curvature of S{r) has the following expansion in powers 
of inverse r: 

TZ = 2r- 2 + TZ- 3 r- 3 + 0(r~ 4 ) , (A4) 

where the coefficient TZ~ 3 may be expressed in terms of the asymptotic shear a of N '. f| 
If we use formula ( |2.1| ) (appropriate for a lightcone embedding) along with the "rest-frame" 
choice 1 = 0, then we have 

k = -\ / 2K. (A5) 

Although this expression differs from ( p.2|) , it is also a lightcone reference, from (|A5| ) we 
find the following radial expansion for the reference mean curvature: 

k\ IBt = -2r- 1 - \r- 2 1l- z + 0(r" 4 ) , (A6) 

which agrees through 0(r~ 2 ) with the k\ IcC expansion obtained in Ref. [Q via the Euclidean 
reference. Hence, all of the results found in Ref. || are also valid for the this choice of 
lightcone reference. In particular, the lightcone-referenced QLE agrees with the TBS energy 
in a suitable null limit. Moreover, in the same limit the "smeared" version of the QLE [which 
incorporates a lapse function into the definition ( J0.1[ )] agrees with Geroch's supermomentum 
(when the latter is evaluated in a Bondi conformal frame). See Ref. [|[ for further details. 
To differentiate between the lightcone and Euclidean references in the large-sphere limit, 
one could examine multipole-moment terms (which arise at higher powers of inverse radius) 
for stationary spacetimes. We hope to return to this issue elsewhere. 



APPENDIX B: CONNECTION AND CURVATURE 



Throughout our discussion, we have taken (— , +, +, +) as the signature of the spacetime 
metric. We have also adopted the index conventions of Ref. |19 for curvature. With these 
choices we shall define the standard Newman-Penrose (np) connection and curvature coeffi- 
cients such that our NP equations match those listed by Newman and Tod in the appendix 
of Ref. |T7| . As we use a somewhat non-standard NP formalism, let us list some of our basic 
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conventions explicitly. Label our null frame as {P, 
name its associated connection coefficients as follows: 



7 
13 
a: 



|( r i2i 
2(^122 
2(^123 

|(ri24 



+ r 43 i) 


K 


= Tl31 


7T = 


1^421 


+ 1^432) 


T 


= Tl32 


V = 


1^422 


+ ^33) 


a 


= Ti33 


P = 


1^423 


+ 1^434) 


P 


= Tl34 


A = 


1^424 



m M } = {ei M , , e 3 M , e 4 M } and 



(Bla,b,c) 

(P,e,f) 
(PTk,li,i) 



(E3,k,l) 



where, for example, r 134 = 1%^^. For the small-sphere limit the null frame is adapted 
to the lightcone N p as described in the first section, and, as a result, the following vanish 
identically for our construction: 



P~ P= P~ P 



K 



r — a — (3 = tt — a — (3 = 0. 



(B2) 



As we retain the freedom to perform r-independent rotations of the space leg m M , we shall 
use those elements of the Geroch-Held-Penrose (ghp) formalism [IS] pertaining both to 
spin- weighted scalars and the "eth" operator 5 on S(r). If a quantity Q transforms as 
Q — > exp(2isx)Q under the rotation m M — > exp(2ix)m fl (with \ independent of r), then Q is 
said to be a spin-weighted scalar of spin weight s (in symbols, SW(Q) = s). As the extents 
of our null normals Z M and n M have been fixed once and for all, we have no need to consider 
the concept of boost weight. 

Next, consider the spacetime Ricci scalar as well as the components of the Weyl and Ricci 
tensors with respect to the null tetrad. With these define the following standard pieces of 
the spacetime curvature: 



C 



4312 J 



(B3a,b,c) 
(BJi,e,f) 
0,h,i) 
®j,k,l) 
&,n,o) 

As mentioned above, with these conventions our NP equations are exactly those listed in the 



#0 


— C*1313 


*1 


= C*1213 


^2 : 


- | (Cl212 


$3 


= Cl242 


^4 


= C*2424 


A = 




$00 


= 2-Rll 


$10 


— 2-^14 


$20 


— 2-^44 


$01 


= 2^13 


$11 


= 4(^12 + -R34) 


$21 


= 2-^24 


$02 


= |-R33 


$12 


= 5-^23 


$22 


= 2-^22 • 



appendix of Ref . |l7j ; however, we shall consider all possible simplification of these equations 
afforded by Eq. ( |B2[ ) and the use of 9 [cf. Eqs. ( |B14[ ) and ( |B16| )]. For the small-sphere limit 
examined in this paper, we consider the pullbacks to the lightcone N p of the curvature 
components ( |B3| ) and assume that each pullback may be expanded as a power series in r 
along the lightcone. That is, we assume 



*i = *° + r^] + r 2 ^ 2 H 

$,, = $°- + r$^- + r 2 $ 2 - + • • • 
A = A + rA 1 + r 2 A 2 + • • • . 



0,1,2,3,4 



i,3 = 0,1,2 



(B4a) 

(He) 



We shall use these expansions for the curvature components along with the radial NP field 
equations in order to obtain radial expansions along N p for those spin coefficients ( |B1| ) used 
in our analysis. We consider the non-vacuum and vacuum cases separately. 
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1. Non- vacuum 



Up to the appropriate order in the affine radius, we have confirmed the non-vacuum 
asymptotic expansions of the spin coefficients given by Kelly et al [0]] and by Dougan. [T^ 
For completeness we recall that this list is 

p = -r- 1 + §r$° + 0(r 2 ) (B4a) 
§r^° + 0(r 2 ) (H|b) 



a 



a = r^a + ±r [a°$° + 2$° - f ° - a°^° ] + 0(r 2 ) 

(3 = -r~ x dP - \r [a°$° 00 - 3^? - a°^° ] + 0(r 2 ) flB|d) 

A = |r [5$° + + 0(r 2 ) (0e) 

A* = -Ir- 1 + \r [tfjj + #° + 2A° + - |$° ] + 0(r 2 ) , (@f) 



where a = y |£. For the scenario at hand the expansions ( |B4| c,d) determine the expansions 

for r and 7r up to 0(r 2 ). Actually, we need only know that A = 0(r) in order to obtain the 
0(r) coefficient for //; however, for completeness we have explicitly given the 0(r) coefficient 
for A. As they are not needed in this paper, we do not list the expansions for 7 and v. 



2. Vacuum 

For the vacuum scenario, we need some of the required spin-coefficient expansions out 
to an order higher than given in either of Refs. |Tl| , [l21| . We obtain the following list: 



p = -r- 1 + ^r 3 ^° + 0(r 4 ) (B5a) 

a = §r^° + \r 2 ^l + |r 3 ^ 2 + 0(r 4 ) flB|b) 
a = r-'a - \r [*° + a ^] - ±r 2 [%\ + a°^ l ] 

+ 3^r 3 [7a°tfgtf{j + 8tf{jtf? - 18a ^ 2 - 35 ^ 2 ] + 0(r 4 ) QB|c) 
/5 = -r~ l oP + ±r [3*S + «°^o] + ^ [ 4 ^i + a °^o] 

-3^r 3 [7a°*°*° - 20tfgtf? - 18« ^ 2 - 158 ^ 2 ] + 0(r 4 ) ([B5|d) 

A = -I r f ° + 0(r 2 ) flB^e) 
H = -Ir" 1 + |r [«§ + *§] + |r 2 [M>^ + 

[sio^^^ - ijSo + 4*?) - p (*°*? + 4* 2 ) - I^ 2 ] + 0(r 4 ) , (Hf) 

where the coefficient 7?. 2 [cf. Eq. (|B21j) 1 found in \i is written out explicitly below in 
Eq. ( |B23| c). Also, the unit-sphere "eth" operator 5 is defined below in Eq. ( |B15[ ). As 



before, Eqs. (|B5|c,d) determine expansions for r and 7r up to 0(r 4 ), and we do not need the 
expansions for 7 and v. 

The particular form of the 0(r 3 ) coefficient p? in the expansion (|B"5|f) for 



p = -|r- 1 + --- + r 3 /i 3 + 0(r 4 ) (B6) 

plays a crucial role in our calculation of the qle's vacuum limit. Therefore, let us sketch 
how to obtain the given form of /i 3 , assuming that we have already determined both p, up 
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to and including 0(r 2 ) and the remaining expansions ( |B5|) as given. Rather than making a 
straightforward appeal to the NP field equations (as we have done to obtain the other spin 
coefficients), we shall instead take advantage of the particular geometry of our construction 
and derive this coefficient via the geometric identity [[18] 



K = -* a + $ u + A - <tA + /i/o . (B7) 

K is the complex Gauss curvature of S(r), and the Ricci scalar of S(r) is simply 1Z = 
2{K + K). Since we work here in vacuo and both p and p are real, Eq. (|B7p implies that 

K = -2(m 2 + ^ 2 )-2(a\ + a\)+App. (B8) 



Now, into this equation insert both Eq. ( |B6|) and the expansions (|B5| a,b,e), and then isolate 



the 0(r 2 ) piece of the resulting expression, in order to establish that 

/** = + + 4*8*8 -J^ 2 - (B9) 

To work this relationship into the desired form, we appeal to the following vacuum NP 
Bianchi identity: 

DV 2 - 3*i = -A^ + tt^i + 3p^ 2 = , (BIO) 

where D = Z^V^ and the action of 5 on \&i is defined as the conjugate of S^i [with 5 as in 
Eq. (|B16| ) and sw(*i) = — 1]. The given form of Eq. flB10| ) is specific to the geometry of 
our construction, but it can be deduced from the (more general) appendix Eq. (A. 4c) given 
in Ref. JIT]. Now, plug the expansions ( pi|) and (|B5|) as well as the radial expansion for 5 



into the Bianchi identity ( |B10| ), thereby obtaining a tower of identities (one identity at each 
power of r). Isolate the particular identity determined at O(r) in the tower, and into this 
equation make repeated substitutions with ( |B24j a,b). These steps lead to 

*! = ^*8*8 + po(*8*? + 4* 2 ), (bii) 



which, upon substitution into Eq. (|B9|), establishes the chosen form of the coefficient p 3 
found in the p expansion (|B5|f). 



3. Intrinsic geometry of S{r) 

In this final appendix subsection we consider the intrinsic geometry of S(r) and collect 
radial expansions along N p for the "eth" operator S, the surface area element dS, and the 
intrinsic curvature scalar 1Z [the Gaussian curvature of S(r) being twice 1Z}. Again, we 
consider the non-vacuum and vacuum cases separately. We obtain the desired intrinsic- 
geometry expansions as follows. First, writing the complex space leg of the null tetrad as 



5 = m^ d/dx" = £ d/d( + r? d/d( , (B12) 

we determine the r expansions along N p for the tetrad coefficients £ and r\ from the NP 
commutator equations, [ff7L|TT| 
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Di = pi + or, (B13a) 
Dr/ = pr/ + a£ (pT3| b) 
5f - = (or - 0)£ - (a - 0)77 , (PI3fc) 

and the radial expansions for the spin coefficients [Eqs. (@a,b,c,d) or Eqs. (p5|a,b,c,d) as 
the case may be]. The given forms of these equations are particular to how our null frame 
has been adapted to the lightcone N p . 

Let us first consider the 9 operator. With both Eqs. (|B"4"|c,d) and the non-vacuum radial 
expansions for £ and rj, we compute the non- vacuum radial expansion for 9, 

d = 5 + s(a-(3) 

= r~% + \r [*°9 - 4s*° + $° 9 + 2s<5° ] + 0(r 2 ) . (B14) 

Here, as above, 5 is the standard NP notation for the directional derivative along m' 1 , s = 
sw(tp) denotes the spin weight of some spin-weighted scalar ip on which 9 acts, and in terms 
of P = 1 + CC we define 

9 = \J\Pd/dQ + 2sa° . (B15) 

A similar calculation shows that the vacuum 9 operator is 

d = 5 + s(a-(3) 
= r~% + |r (\I/°6 - 4s*?) + ^r 2 (*J8 - 5s**) 

+ 3^r 3 [7*°*°9 + 18* 2 9 - 18s (9 * 2 ) - 125*°*°] + 0(r 4 ) , (B16) 

with 5, s, and 9 as before. 

Next, we substitute the radial expansions we find for £ and rj into area element of S(r), 



dS = UQP 2 ( 



(B17) 



Here dQ is the standard unit-radius round sphere area element. For the non- vacuum scenario 

dS = dttr 2 [1 - §r 2 $° + 0(r 3 )] , (B18) 

while for the vacuum scenario 

dS = dttr 2 [1 - ^r 4 *°*° + 0(r 5 )] . (B19) 

These expansions may be checked with the well-known formula D(dS) = —2pdS for the 
change in the surface- area element. fllQfl 

Let us turn to the Gauss curvature of S(r). An order-by-order examination of the identity 

(99 - 99V = \K<P (B20) 

[where Sw(</?) = 1, but ip is otherwise arbitrary, and we assume the expansion <p = ip° + 
np 1 + r 2 ip 2 + ■ ■ •] determines the following expansion for the intrinsic curvature of S(r): 
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K = 2r~ 2 + 11° + rTZ 1 + r 2 Tl 2 + 0(r 3 ) . (B21) 



In examining Eq. (|B20 ), we must use the appropriate expansion for 9, Eq. (|B14j ) or Eq. ( |B16| 



as the case may be. For the non-vacuum scenario, we find the first relevant coefficient to be 

K° = | [$° + g $? + e $? - 29 f J - 29 *°] . (B22) 

For the vacuum case we need more coefficients. Tedious but straightforward calculation 
yields the set: 

TZ° = -§ (5 *° + Bo*?) (B23a) 

n 1 = -| (s *i + s *y @b) 

n 2 = - \m - \m - gs (*>?) - p (*g*») . 



We remark that, along the way to obtaining these coefficients of 1Z from Eq. ( |B20| ), one finds 
terms involving derivatives of the coefficients of ip. However, one knows that such terms must 
vanish, as the right-hand side of Eq. ( |B20[ ) involves no derivatives of <p. Therefore, one may 



simply discard these terms. A detailed examination (order by order in r) of the Bianchi 
identities (non- vacuum or vacuum as the case may be) shows that such terms indeed vanish. 



4. Vacuum Bianchi Identities 

Consider the vacuum Bianchi identities as given in in the appendix Ref. [pjj. For the 



situation at hand, we may simplify these equations by writing them in terms of S rather 

than 5, and we do so. Then, considering the equations order by order in r, we find the 
following at the lowest order in r: 

5 ^o = 8 *° = 4*? (B24a,b) 

g o $o = _i*g g o *° = 3*° (jB2gs,d) 

s *° = -*° 8 *° = 2*° B,f) 

g *° = -§*° 8 *° = *° (iHg,h) 

g *° = -2*o e *° = o. (gg§j) 

from these one can derive results such as S 8 *2 — ~ 3*°- 
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FIGURES 





FIG. 1. Geometry of the Lightcone Reference. 
Respectively, the top and bottom boxes depict the physical and reference spacetimes. In the top 
box the shaded surface spanning S(r) is the 3-surface S, while in the bottom box the shaded surface 
spanning S(r) is the 3-surface Z. Whether viewed as the intersection N p f)Y] or the intersection 
N g P|£, the 2-surface S(r) has the same intrinsic 2-metric. Our limit construction gives us X, 
but we must choose Z; and, moreover, our choice of Y. must be determined solely by the intrinsic 
2-metric of S(r). Our choice and its physical motivation are described in Subsection 2. A. However, 
we note here that, whenever S(r) is at all distorted from perfect roundness (as it generally will 
be) , Y. is not flat Euclidean 3-space E 3 (because the intersection of E 3 with the genuine lightcone 
N 9 would be a round sphere). Choosing such a lightcone reference, we assign the zero value of 
the energy to that (shaded) portion of Y. contained within S(r), and compute the energy of (the 
shaded portion of) £ relative to this zero-point. 
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